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ABSTRACT 



> 

' We present results of large N— body-hydrodynamic simulations of galaxy forma- 

ly-j , tion. Our simulations follow the formation of galaxies in cubic volumes of side 100 Mpc, 

■ in two versions of the cold dark matter (CDM) cosmogony: the standard, Q = 1 SCDM 

model and the flat, f2 = 0.3 ACDM model. Over 2000 galaxies form in each of these 

, simulations. We examine the rate at which gas cools and condenses into dark matter 

■ halos. This roughly tracks the cosmic star formation rate inferred from observations at 

various redshifts. Galaxies in the simulations form gradually over time in the hierar- 

O ■ chical fashion characteristic of the CDM cosmogony. In the ACDM model, substantial 

I 1 galaxies first appear at z ~ 5 and the population builds up rapidly until z = 1 after 

which the rate of galaxy formation declines as cold gas is consumed and the cool- 
ing time of hot gas increases. In the SCDM simulation, the evolution is qualitatively 
similar, but it is shifted towards lower redshift. In both cosmologies, the present- 
day K-band luminosity function of the simulated galaxies resembles observations. The 
galaxy autocorrelation functions differ significantly from those of the dark matter. At 
the present epoch there is little bias in either model between galaxies and dark matter 
on large scales, but a significant anti-bias on scales of ~ 1 h^ 1 Mpc and a positive 
bias on scales of ~ 100/i _1 kpc . The galaxy correlation function evolves little with 
redshift in the range z = — 3, and depends on the luminosity of the galaxy sample. 
The projected pairwise velocity dispersion of the galaxies is much lower than that of 
the dark matter on scales less than 2 h~ x Mpc . Applying a virial mass estimator to the 
largest galaxy clusters recovers the cluster virial masses in an unbiased way. Although 
our simulations are affected by numerical limitations, they illustrate the power of this 
approach for studying the formation of the galaxy population. 

Key words: cosmology: theory — galaxies: formation — galaxies: kinematics and 
dynamics — hydrodynamics — methods: numerical 



1 INTRODUCTION 

A detailed understanding of galaxy formation is one of the 
central goals of contemporary astrophysics. Over the past 
decade, this goal seems, finally, to have come within reach. 
On the observational side, data from the Keck and Hubble 
Space telescopes have revolutionised our view of the high 
redshift Universe. From the theoretical point of view, it is 



clear that studying formation involves a synthesis of ideas 
from a wide range of specialities. A full theoretical treat- 
ment requires consideration of the early Universe processes 
that created primordial density fluctuations, the non-linear 
dynamics that result in the formation of dark matter halos, 
the dissipational processes of cooling gas, the microphysics 
and chemistry that precipitate star formation, the feedback 
due to energy exchange between supernovae and interstellar 
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gas, the merging of galaxy fragments and the interaction of 
galaxies with their large-scale environment. 

Because of its non-linear character, lack of symmetry 
and general complexity, galaxy formation is best approached 
theoretically using numerical simulations. As a minimum, 
realistic simulations must be able to follow the evolution of 
the dark matter from appropriate initial conditions, together 
with the coupled evolution of a dissipative gas component 
that will eventually produce the visible regions of galaxies. 
The advent of large computers, particularly parallel super- 
computers, together with the development of efficient algo- 
rithms, has given considerable impetus to this approach in 
recent years. 

Both Eulerian and Lagrangian gas dynamics techniques 
have been used to simulate galaxy formation. Currently, only 
the latter have achieved sufficient dynamic range to resolve 
galaxy formation within cosmological volumes. For example, 
the recent large Eulerian simulations of Blanton et al. (1999) 
have spatial resolution of order 300 — 500 kpc, whereas typ- 
ical Lagrangian techniques can readily achieve spatial reso- 
lutions of order 20 kpc or better (e.g. Carlberg, Couchman 
& Thomas 1990, Katz, Hernquist & Weinberg 1992, 1999, 
Evrard, Summers & Davis 1994, Frenk et al. 1996, Navarro 
& Steinmetz 1997). Recently, the power of the traditional, 
fixed-grid Eulerian approach has been greatly enhanced by 
the inclusion of adaptive mesh refinement techniques which 
have been applied to the study of the early phases of galaxy 
formation (Abel et al. 1998). 

In this paper, we carry out simulations of galaxy for- 
mation using the Lagrangian technique of "Smooth particle 
hydrodynamics" or SPH (Gingold & Monaghan 1977; Lucy 
1977). We focus on the thermodynamic history of the gas 
that cools to make galaxies, on the rate at which galaxies 
build up with time, on their abundance at the present day 
and on their spatial and velocity distributions. Early results 
from one of the simulations that we analyse here were pre- 
sented in a previous paper (Pearce et al. 1999), where we 
focused on the spatial distribution of galaxies. Our work fol- 
lows on from earlier attempts to grow virtual galaxies in a 
cosmological volume using the SPH technique, most notably 
by Carlberg, Couchman & Thomas (1990), Katz et al. (1992, 
1996, 1999), Evrard et al. (1994) and Frenk et al. (1996). 
These simulations succeeded in resolving individual galaxies 
and allowed the first investigations of their distribution and 
of the dynamics of galaxies in clusters. However, the vol- 
umes modelled in these early simulations were rather small, 
permitting only limited statistical studies. Our simulations 
contain about one order of magnitude more particles and 
produce about 40 times more galaxies than the largest pre- 
vious simulation of this kind by Katz, Weinberg & Hernquist 
(1996). SPH simulations with even higher resolution have 
been used to study the formation of individual galaxies (e.g. 
Navarro & Steinmetz 1997, Weil, Eke & Efstathiou 1998). 

Simulations of the kind studied here are complemen- 
tary to semi-analytic studies of galaxy formation (e.g. Kauff- 
mann et al. 1997, 1999a,b; Diaferio et al. 1999; Guiderdoni et 
al. 1998; Benson et al. 1999a,b; Somerville & Primack 1999, 
Cole et al. 2000.) The main advantage of the hydrodynamic 
simulations is that they follow the dynamics of diffuse cool- 
ing gas in full generality whereas the semi-analytic treat- 
ment involves a variety of approximations such as spherical 
symmetry and quasistatic evolution. A comparison of re- 



sults from the two methods has recently been performed 
by Benson et al. (2000c). Both approaches require a phe- 
nomenological treatment of critical processes such as star 
formation and feedback since these involve scales far below 
the resolution limit of all current simulations. 

The remainder of this paper is laid out as follows. In 
Section 2 we describe our simulations and the physical ap- 
proximations we have made. In section 3 we discuss the evo- 
lution of the global properties of the gas, the extraction of 
a dark matter halo catalogue, and the properties of the gas 
within these halos. In Section 4 we focus on the properties 
of galaxies, their formation histories, luminosity function, 
age distribution, etc. In Section 5 we consider the correla- 
tion function of the galaxies as a function of redshift and 
galaxy mass, the relative pairwise velocity distributions of 
galaxies and dark matter, and we determine how well the 
galaxies trace the mass distribution within the largest dark 
matter halos. We conclude in Section 6. This work is part 
of the programme of the Virgo consortium for cosmological 
simulations. 



2 SIMULATIONS 

Our simulations were performed using a parallel, adaptive, 
particle-particle, particle-mesh SPH code (Pearce & Couch- 
man 1997) which can be run in parallel on the Crays T3D 
and T3E, or serially on a single processor workstation. This 
code is essentially identical in operation to the publicly re- 
leased version of hydra, described in detail by Couchman, 
Thomas & Pearce (1995). 

We have simulated galaxy formation in two cold dark 
matter (CDM) models (e.g. Davis et al. 1985) with the same 
values of the cosmological parameters as the SCDM and 
ACDM models of Jenkins et al. (1998). For SCDM these are 
as follows: mean mass density parameter, Qq = 1; cosmolog- 
ical constant, A/(3Hq) — 0; Hubble constant (in units of 100 
km s -1 Mpc -1 ), h = 0.5; and rms linear fluctuation ampli- 
tude in 8/i -1 Mpc spheres, as = 0.6. For ACDM , the cor- 
responding parameter values are: Qo = 0.3; A/{3H$) = 0.7; 
h = 0.7; and a 8 = 0.9. Each simulation followed 2097152 
dark matter particles and 2097152 gas particles in a cube 
of side 100 Mpc and required around 10000 timesteps (and 
~ 10 processor hours on a Cray-T3D) to evolve from z = 50 
to z = 0. In both cases, we employed a /3-spline gravi- 
tational softening which remained fixed at the Plummer- 
equivalent comoving value of 50 kpc until z = 1.5 (z = 2.5 
for ACDM ). Thereafter, the softening remained fixed at 
50/(1 + z) — 10ft -1 kpc in physical coordinates, and the 
SPH smoothing length was set to match this value. 

We deliberately chose to set the gas mass per particle 
to be ~ 2 x 10 9 Mq in both cosmologies so as to have iden- 
tical mass resolution. As we typically smooth over 32 SPH 
neighbours, the smallest resolved objects have a gas mass of 
- 6.4 x 10 10 M Q . The baryon fraction, Q b , was set from the 
nucleosynthesis constraint, Qbh 2 = 0.015 (Copi, Schramm 
& Turner 1995) and although this is somewhat lower than 
the more recent limits given by Tytler et al. (2000), it was 
the current value at the time these simulations were carried 
out. We also assumed a constant gas metallicity of 0.3 times 
the solar value in both cases. We expect more gas to cool 
globally in the ACDM model because structure forms earlier 
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Figure 1. An example of excessive cooling near a density interface. The two panels on the left show radial temperature profiles and the 
two panels on the right radial density profiles for a galaxy surrounded by a hot gaseous corona. In both cases, the top panel shows the 
result of decoupling the cold and hot gas phases as described in the text, while the bottom panels shows the profiles with no decoupling 
of the cold and hot phases. In the latter case, the density associated with hot particles near the interface between the two phases is 
overestimated, leading to enhanced cooling. The vertical dashed line shows the mean search length for particles which have cooled into 
the galaxy and the vertical dotted line the mean search length for the hot halo particles. In the case in which phase decoupling has been 
implemented, the particles of the central galaxy are shown as dots rather than triangles. 



in this model. With our chosen parameters, our simulations 
were able to follow the cooling of gas into galactic-sized dark 
matter halos. The resulting "galaxies" are typically made 
up of 50 — 1000 particles and so have cold gas masses in 
the range 10 11 M - 2 x 10 13 M . With a spatial resolution 
of 10/i _1 kpc, we cannot resolve the internal structure of 
these galaxies and we must be cautious about the possibil- 
ity of enhanced tidal disruption, drag, and merging within 
the largest clusters. However, as we argue below, there is no 
evidence that this is a major problem. 

As in all numerical simulation work, approximations 
and compromises are necessary if the calculation is to com- 
pleted within a reasonable amount of time. An important 
choice is the effective resolution of the simulation. Increasing 
the spatial resolution requires increasing both the number 
of particles (to prevent 2-body effects) and the number of 
timesteps (to follow smaller structures). Our chosen value of 
10/i _1 kpc for the gravitational softening is larger than the 
scalelength of typical galaxies, but is the best that we can 
do with our current computer resources. For a given simu- 
lation volume, the number of particles determines the mass 
resolution. In cosmological simulations that involve only col- 
lisionless dark matter, the choice of particle number is often 
driven by computer memory limitations. Cosmological gas 
dynamics simulations are more computationally demanding 
than collisionless simulations and memory considerations do 
not, at present, limit their size. Since approximately 10% of 
all the baryons in the universe have cooled by the present to 
form stars and dense galactic gas clouds (Fukugita, Hogan 
& Peebles 1998), and around 50 particles are required to 
resolve a galaxy in a simulation, a simple estimate of the 



number of gas particles required to obtain N ga i galaxies is: 
10 x Ng a i x 50, or roughly half a million particles for 1000 
galaxies. In practice, many galaxies contain more than the 
minimum number of particles and so this is an underesti- 
mate. Since galaxies of the characteristic luminosity, L* or 
greater, have a space density of ~ 0.01/(/i -1 Mpc) 3 , a simu- 
lation of a volume like those studied by Jenkins et a!. (1998), 
of side ~ 250/i _1 Mpc, would produce over 100, 000L* galax- 
ies and require around 10 s gas particles. For this work, we 
have chosen to use just over 2 x 10 6 million gas particles in 
a volume of side only 100 Mpc. 

In addition to the choice of simulation parameters, other 
choices need to be made to model physical processes that 
operate on sub-resolution scales. Star formation and the as- 
sociated feedback from energy released in supernovae and 
stellar winds are the most obvious sub-resolution processes 
relevant to simulations of galaxy formation. As the reso- 
lution is increased, progressively smaller, denser halos are 
resolved whose gas has a short cooling time. In hierarchi- 
cal clustering theories, a simple model predicts that most of 
the baryonic material should have cooled (and presumably 
turned into stars), in small, subgalactic objects at high red- 
shift (White & Rees 1978, Cole 1991, White & Frenk 1991). 
Since this does not appear to be the case in the universe, 
some mechanism must have prevented much of the gas from 
cooling. Reheating of the gas by the energy released in the 
course of stellar evolution is a commonly invoked mechanism 
to quench star formation in small galaxies. This, however, 
is still a very poorly understood process. Various simplified 
models have been implemented in simulations of galaxy for- 
mation (e.g. Katz 1992; Navarro & White 1993; Metzler & 
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Evrard 1994; Steinmetz & Muller 1995; Evrard, Metzler & 
Navarro 1996; Katz et al. 1996; Mihos & Hernquist 1994, 
1996; Gerritsen & Icke 1997; Navarro & Steinmetz 1997). 

In this paper, we have taken the simpler route of ne- 
glecting feedback effects altogether. Instead, we appeal to 
the resolution threshold of the simulation to prevent all the 
gas from cooling. This is analogous to assuming a model in 
which gas cooling is suppressed within dark matter halos of 
mass below ~ 6.4 x 10 10 Mq), whilst leaving larger halos 
unaffected. Although relatively simple, this model gives re- 
sults that are not too dissimilar from those produced by de- 
tailed semi-analytic models in which feedback is treated with 
greater care. In semi-analytic modelling, the star formation 
and feedback prescriptions are tuned by the requirement 
that the model should reproduce observable properties of the 
local galaxy population such as the luminosity function. The 
difference in the amounts of gas that cool in our simulations 
and in the semi- analytic model of Cole et al. (2000) are illus- 
trated in Fig. 1 of Pearce et al. (1999). Near the resolution 
threshold, approximately 50% more gas cools in the simu- 
lation than in the semi-analytic model and the difference 
falls rapidly with increasing galaxy mass. A detailed com- 
parison of results from these simulations and semi-analytic 
modelling may be found in Benson et al. (2000c). 

The second area in which sub-resolution physics is im- 
portant is the interaction of cool gas clumps with the sur- 
rounding hot gas. In reality, the gas within such clumps 
would presumably become incorporated into molecular 
clouds or turned into stars. These processes occur on scales 
that are several orders of magnitude below the resolution 
limit of current simulations and so must be modelled phe- 
nomenologically. Different authors have chosen to do this in 
different ways. For example, Evrard et al. (1994) and Frenk 
et al. (1996) opted to leave the cold gas alone. Although this 
is the simplest possible approach, it has the disadvantage 
that the smoothing inherent in the SPH method can arti- 
ficially boost the cooling rate of hot gas that comes into 
contact with the gas that has cooled to very high density in- 
side a galaxy. This process is illustrated in the lower panels 
of Fig. 1 which show the temperature and density profile of 
a galaxy in the simulation. Model galaxies are typically sur- 
rounded by hot gaseous coronae. At the interface with the 
cold galactic gas, the density of the hot coronal gas is over- 
estimated and its cooling rate is artificially enhanced. This 
gives rise to the rapidly cooling material which is clearly vis- 
ible in the region between the effective search lengths of the 
cold and hot material, indicated by the dashed and dotted 
lines in the figure respectively. The artificial cooling at this 
contact interface arises because SPH is not designed to deal 
with the steep density gradients that develop in a multiphase 
medium. As Thacker et al. (2000) have shown, this limitation 
is present to some extent in all standard implementations of 
SPH, although recently Ritchie & Thomas (2000) have pro- 
posed an alternative formulation of SPH utilising the gas 
pressure rather than the gas density which overcomes this 
limitation. 

An alternative approach to leaving the galactic cold 
dense gas alone is to replace it with collisionless "star" par- 
ticles, thus decoupling it entirely from the rest of the bary- 
onic material (Navarro & White 1993, Steinmetz & Muller 
1995, Katz et al. 1996). While this might appear as a more 
realistic solution than leaving cold material as gas, in prac- 
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Figure 2. The effect of decoupling the hot halo gas from the 
cold galactic material. The graph shows the ratio of the masses of 
matched clumps found in two test simulations with similar param- 
eters and resolution to those assumed in our main simulations. In 
one of the test simulations, the hot gas was decoupled from the 
cold gas, as described in the text; in the other it was not. The 
quantity plotted along the j/-axis is the ratio of the clump mass 
in the first to the clump mass in the second of these test simula- 
tions. Only the largest object (which contains over 1000 particles 
in the simulation without decoupling) is significantly affected. 



tice, the relatively large softening lengths normally used in 
SPH simulations imply unrealistically low binding energies 
for the star particles. As a result, stellar objects are prone 
to disruption by collisions and tidal encounters, particularly 
in strongly clustered regions. 

In this study, we have chosen an approach which is in- 
termediate between the two extremes just discussed. In our 
phenomenological model, hot gas, defined to be hotter than 
10 5 K, is prevented from interacting with (galactic) material 
at a temperature below 12000 K (but not vice- versa). All 
other SPH forces remain unchanged. The threshold tem- 
perature of 10 K is less than the virial temperature of the 
smallest resolved halos in our simulations. In our implemen- 
tation of SPH, this assumption preserves force symmetry 
and momentum conservation, and effectively prevents over- 
cooling of hot gas induced purely by the presence of a nearby 
clump of cold, dense gas. This is clearly seen in the top pan- 
els of Fig. 1. In contrast to the standard model in which the 
hot and cold phases remain strongly coupled, there is now 
no interface discontinuity. As a result, the density profile in 
the vicinity of a cold clump remains smooth and a "cooling 
flow" develops around it. Recently this effect has been con- 
firmed by Croft et al. (2000) who suggest several alternative 
routes to circumvent it. 

In our scheme, the cold gas "feels" all the usual SPH 
forces, including viscous drag when it moves in a hot inter- 
galactic medium, and interacts with other cold gas clouds in 
the usual way. The outer accretion shock of a galaxy cluster 
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Figure 3. The relative fractions of gas in each of three phases: 
(i) cold, galactic gas (temperature less than 12000 K, overden- 
sity greater than 10); (ii) hot halo gas (temperature greater than 
10 5 K); (iii) cold, diffuse gas in "voids" (overdensity less than 10) . 
The bold lines correspond to the SCDM simulation and the thin 
lines to the ACDM simulation. 



Figure 4. The rate at which gas cools, per unit volume, in the 
SCDM simulation as a function of time since the Big Bang. The 
data points are inferred star formation rates in the real universe, 
from Steidel et al. (1999). The match between the simulation re- 
sults and the data indicates that, in spite of resolution effects 
and other limitations, the simulations provide a reasonably real- 
istic model. 



is still properly modelled as the infalling material can "see" 
the hot gas in the cluster. Our approximation has only a 
small effect on the masses of small galactic objects as may 
be seen in Fig. 2, where we compare the masses of galax- 
ies formed in two small test simulations. These had similar 
parameters and mass resolution to those employed in our 
main simulations, but in one case the hot gas was decou- 
pled from the cold gas and in the other case it was not. For 
most galactic-sized objects, this change produced a slight 
reduction in mass, but for the largest object, which has over 
1000 particles without decoupling, the mass is considerably 
reduced to 370 particles. This effect, as well as the results of 
varying several simulation parameters are explored further 
by Kay et al. (2000). 



3 THE PROPERTIES OF THE GAS 

In this section, we examine both the evolution of the gas in 
the simulation as a whole and the properties of gas within 
virialised dark matter halos. 

3.1 Global properties of the gas 

Following Kay et al. (2000), we define three different gas 
phases: (i) cold galactic gas (temperature below 12000 K and 
overdensity greater than 10) ; (ii) hot halo gas (temperature 
above 10 5 K) ; and (iii) cold, uncollapsed gas (overdensity 
less than 10; this is essentially diffuse gas in "voids.") We 
ignore the effects of photoionization because at this mass 
resolution a photoionising field has a negligible effect upon 
galactic or halo properties. In practice, little gas is both cold 



and at overdensities between 5 and 1000, so that effectively 
the uncollapsed and galaxy phases are completely disjoint. 
The evolution of the fraction of gas in each of these phases is 
illustrated in Fig. 3 for the two cosmological models that we 
have calculated, SCDM (bold lines) and ACDM (thin lines). 
These gas fractions are very similar to those derived using 
a semi-analytic model of cooling gas implemented in these 
simulations by Benson et al. (2000c). 

The properties of cold, uncollapsed material in voids are 
only followed in an approximate manner in our simulations. 
The SPH method generally has poor spatial resolution in re- 
gions in which the gas is diffuse because it must smooth over 
a sparse particle distribution. This limitation is exacerbated 
in our specific implementation which imposes a maximum 
search length, resulting in a minimum resolvable density of 
a few times the mean density in the computational volume. 
Test simulations without such a minimum resolvable den- 
sity demonstrate that these inaccuracies do not affect the 
properties of the galaxies because the gas forces are small 
in void regions and have no effect on the dense regions that 
ultimately form galaxies (Kay et al. 2000). However, out- 
wardly propagating shocks generated by the formation of 
large clusters are poorly followed in underdense regions and 
so we can say little about the detailed properties of voids 
(beyond their mere presence). This is a regime in which Eu- 
lerian techniques are more accurate than SPH {e.g. Cen & 
Ostriker 1996). 

As the evolution proceeds the cold, diffuse phase is de- 
pleted as gas falls into dark halos where it is shocked into 
the hot phase, thereafter cooling into galaxies. The time 
derivative of the dashed curve that tracks the galaxy phase 
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Figure 5. A comparison of the cumulative mass function of dark 
halos in the ACDM gas + dark matter simulations and an other- 
wise identical pure dark matter simulation. The halo masses of the 
gas + dark matter simulation have been scaled by f2o/(^0 — Qb)- 
For this plot the halos were identified using a friends-of-friends 
algorithm with a linking length of 0.2. The shaded region shows 
the expected number of halos, allowing for lcr Poisson errors be- 
cause of the finite simulation volume, as predicted by the formula 
given in Jenkins et al. (2000). 



in Fig. 3 effectively defines the global star formation rate 
in the computational volume. This is plotted in Fig. 4. The 
star formation rate rises rapidly at early times, has a broad 
maximum around z ~ 1 — 2 and then gently declines to the 
present. This behaviour is broadly similar to that seen in 
real data (Steidel et al. 1999), indicating that the limited 
resolution of the simulations and other approximations do 
not result in too unrealistic a model. At late times, Fig. 3 
shows that there is a substantial transfer of gas from the 
uncollapsed phase into the hot phase, with around 50% of 
the gas ending up at temperatures above 10 5 K in both cos- 
mological models. At z ~ 1, when only around 30% of the 
gas is above 10 5 K, the gas fractions in the different phases 
in our simulations agree well with those in the simulations 
of Evrard et al. (1994) (which were stopped at this point.) 

As expected, more gas cools into the galaxy phase in 
the ACDM than in the SCDM model. This is due to our 
adoption of a single gas particle mass in the two simulations 
which results in the same effective resolution in both cases. 
More gas cools at early times in the ACDM simulation be- 
cause structure forms earlier in this case. If all the gas that 
cools into galaxies were assumed to turn into stars, then the 
SCDM model would have roughly the observed mean mass 
density in stars. In the ACDM model, on the other hand, 
about twice as much gas cools and, if all this gas were as- 
sumed to turn into stars, an uncomfortably large value of 
the mass-to-light ratio would then be required to match the 
observed mean stellar mass density. The amount of gas that 
cools in our simulations is determined by the mass resolution 



and a small decrease in the resolution of the ACDM model 
would reduce the final amount of cold gas to the same level 
as in the SCDM model. 

3.2 Properties of gas in halos 

To construct a catalogue of dark matter halos from the sim- 
ulations, we first identified cluster centres using the friends- 
of-friends grouping algorithm (Davis et al. 1985) with a small 
value of the linking length, b = 0.05. We then grew a sphere 
around the centre of mass of each of these halos until the 
mean overdensity within it reached a value of 178 in the 
case of SCDM or 324 in the case of ACDM . These are the 
overdensities corresponding to the virial radius according to 
the spherical top-hat model for clusters (Eke, Cole & Frenk 
1996.) The resulting set of halos was cleaned by removing 
the smaller of any two overlapping halos. This procedure 
produced catalogues of 1320 halos in SCDM and 1795 ha- 
los in ACDM , with more than 28 particles per halo in each 
case. 

To check whether the presence of baryons affects the 
halo catalogue, we ran a pure N-body simulation of the 
ACDM model and extracted a halo catalogue from both 
it and the baryonic simulation with gas cooling. The mass 
function of halos in the test simulation is almost indistin- 
guishable from that in the original simulation, as illustrated 
in Fig. 5. This figure also compares the mass functions in 
the simulations with the mass function given in Jenkins et 
al. (2000). This gives a better fit than the predictions of the 
Press-Schechter model (Press & Schechter 1974). Although 
for the values of Qb that we have assumed here, the pres- 
ence of gas has little effect on the halo mass function, it 
does affect the central structure of the halo into which it 
cools (Pearce et al. 2000). 

In Fig. 6 we plot the mass fraction of gas within the 
virial radius of the halos normalized to the mass fraction in 
the simulation as a whole, together with the proportions in 
the hot and cold phases. The total gas mass fraction varies 
little with halo mass over the range 1O 12 M <M<1O 15 M , 
although there is some indication of an upturn at low masses 
in the ACDM case and at high masses in the SCDM case. 
The mean value of the gas fraction is ~ 0.85 for ACDM and 
~ 0.8 for SCDM with a scatter of ~ 0.05, indicating that 
the gas is slightly more extended than the dark matter in 
both models. A few small mass halos have a gas mass frac- 
tion greater than unity. These halos, however, are near the 
resolution limit of the simulations. The crosses show the 
mass fraction of hot (T > 12000 K) gas in the halos. This 
increases rapidly with halo mass, from ~ 10 — 20% for the 
smallest halos in the simulations to ~ 60 — 80% for the 
largest. The fraction of cold (T < 12000 K) gas shows the 
opposite trend and varies from ~ 70 — 80% at the small mass 
end to ~ 10 — 20% at the high mass end. 

The fraction of baryonic material that cools within dark 
matter halos can also be calculated using semi-analytical 
techniques (White & Frenk 1991). As we discussed in the 
preceeding section, in semi-analytic models the cooling of 
gas is regulated by feedback processes due to energy liber- 
ated by stars. In these simulations we ignore feedback al- 
together and cooling in small mass objects is then limited 
by resolution effects. This leads to some differences in the 
gas fractions derived with these two techniques. For exam- 
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Figure 6. The ratio of gas to dark matter mass (xQ c( j m /f2(,, where £l c dm denotes the dark matter contribution to Qq) within the virial 
radius of ~ 1500 dark matter halos found in the simulations. The circles give the total gas fraction, the triangles the fraction of gas 
below 12000 K and the crosses the fraction of gas above 12000 K. The solid line indicates the effective resolution limit of 32 gas particles 
within a halo. Panel (a) corresponds to the ACDM simulation and panel (b) to the SCDM simulation. 
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Figure 7. The relation between the mean, mass-weighted temperature of hot (> 10 5 K) gas and the mass of its host dark matter halo. 
The dashed line shows the theoretical expectation, M oc T5 , of a simple equilibrium model for the halos. The tight correlation between 
T and M indicates that the gas is close to equilibrium in all but the smallest halos in the simulations (where the temperatures can be 
overestimated due to the proximity of a larger halo, see text). 



pie, in the semi-analytic model of Cole el al. (2000), the gas 
fraction of cold gas has a maximum of ~ 40% at a halo 
mass of 1O 12 M0. This is close to the resolution limit of 
our simulations. At small masses, the cold fraction is sup- 
pressed by feedback and at larger masses, it is suppressed 



by the long cooling time of the diffuse gas in halos. The first 
of these effects is not present in our simulations (which do 
not resolve the regime where the feedback is important in 
the semi-analytic models), but the second one is. On these 
large mass scales, the simulations agree well with the semi- 
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Figure 8. Evolution of the abundance of "galaxies" in the simulations as a function of their cold gas mass. The mass functions are 
plotted at the redshifts given in the legend. The left panel shows results for ACDM and the right panel for SCDM. 



analytic results and even over the entire range of mass, the 
difference is only of the order of 50% (Benson et al. 2000c). 
Overall, somewhat more mass cools in the simulations than 
in the semi-analytic model, suggesting that feedback effects 
are not negligible even above the resolution limit of our sim- 
ulations. 

The hot gas that falls into a halo is quickly thermalized 
by shocks and heated to the virial temperature. This is illus- 
trated in Fig. 7 which shows the relation between the mean, 
mass-weighted gas temperature (for gas with T > 10 5 K) 
and the virial mass of the host halo. For masses above a few 
times 10 12 Mq, there is a very tight correlation between gas 
temperature and halo mass, indicating that the gas is re- 
laxed and close to equilibrium. At the small mass end, there 
are some halos with anomalously large temperatures. These 
tend to be in the vicinity of larger halos whose own hot 
halos contaminate the smaller ones, boosting their temper- 

3 

ature. The dashed line is the relation M ocT? predicted by 
a simple, spherically symmetric, equilibrium model for the 
halo and its gas. The simulations follow this relation well 
over most of the mass range plotted in the figure. The X-ray 
properties of the largest halos in our simulations have been 
studied by Pearce et al. (2000). 



4 THE PROPERTIES OF GALAXIES 

We now consider the properties of the population of galax- 
ies in the simulations, focussing on their abundance, ages, 
merger histories, and individual star formation rates. 

4.1 Demographics of the galaxy population 

We identify "galaxies" in the simulations with dense clumps 
of cold gas. Locating them is straightforward except in the 



small number of cases when they are undergoing a merger 
or are being tidally disrupted within a large cluster. To find 
galaxies, we use the same friends-of-friends algorithm that 
we employed for finding dark halos, but with a linking length 
of only b = 0.02 for SCDM and b = 0.0134 for ACDM, values 
that are 10% of those required to obtain virialised halos in 
these cosmologies. When calculating the mass of each galaxy, 
we consider only particles with temperature below 12000 K, 
a condition that rejects only very few particles. In practice, 
as the cold clumps of gas that make up the galaxies have 
contracted by a large factor the galaxy catalogues and the 
masses and positions of each are insensitive to these choices. 
We find a total of about 2000 resolved galaxies in each of 
the two simulations at the present day. 

The evolution in the abundance of galaxies is shown in 
Fig 8. Here, we characterize each galaxy by its cold gas mass 
and plot the cumulative abundance at various epochs in our 
two simulations. In the ACDM model, the first substantial 
galaxies begin to form at a redshift z ~ 5, although some 
fairly massive objects are present even before this. The mas- 
sive galaxy population builds up rapidly between z ~ 5 and 
z ~ 3 when it increases approximately ten- fold. A similar 
increase occurs between z ~ 3 and z ~ 1. As time pro- 
gresses, the mass of the largest galaxies also increases. After 
z = 0.5, the growth in the population slows down somewhat 
because the rate at which new galaxies are forming is more 
or less balanced by the rate at which existing galaxies are 
destroyed by mergers and tidal disruption. At z = 0, the 
abundance of galaxies with gas mass M > 10 11 Mq approx- 
imately matches the observed abundance of L* or brighter 
galaxies of ~ 0.003 per /i~ 3 Mpc 3 . The evolution in the 
SCDM model is qualitatively similar, but shifted towards 
lower redshifts. 

The creation sites of galaxies at recent times are illus- 
trated in Fig. 9 which shows the locations where new gas 
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Figure 9. The circles indicate the locations where new gas cools 
after z = 0.3 in the SCDM model, scaled in radius by the mass 
of new cold material. The grey shading indicates the underlying 
dark matter density at z = 0. 



cools between z = 0.3 and the present in the SCDM model. 
(Broadly similar behaviour is seen in the ACDM model.) 
The radius of the circle around each of these locations is 
scaled by the mass of new cold material. Some gas cools 
onto preexisting objects (note the small circles at the centre 
of some clusters), but 383 new resolved galaxies, about 15% 
of the total, formed during this period. The new galaxies 
are small and tend to form along filaments or in low density 
regions, away from the larger clusters. 

Galaxies are destroyed preferentially in large clusters, 
where tidal forces are largest. For example, about 20% of the 
cold material within the virial radius of the largest halo at 
z — 0.3 is reheated before z — 0, and is associated with the 
destruction of about half the small (around 1/3L*) galaxies 
in this halo. As discussed earlier, tidal disruption of galaxies 
in clusters is likely to be enhanced by the limited spatial 
resolution of the simulations which results in artificially low 
binding energies for small objects. Larger simulations are 
required to investigate this important problem in detail. 

The major consequence of "decoupling" the hot gas 
from the cold phase, as discussed in the preceeding section, 
is that no supermassive objects form. The biggest galaxy 
in the ACDM model contains 5 x 10 12 M of cold gas. The 
turnover in the abundance of galaxies at small masses in the 
two simulations is a numerical artifact. At least 32 particles 
(corresponding to a mass of 6.4 x 10 10 M Q ) of gas must be 
present for cooling to be efficient; objects below this are in 
the process of forming or are being ablated as they move 
in a hot halo. The effective mass threshold in both simula- 
tions is actually slightly higher than the nominal value and 
corresponds to a mass threshold around 10 11 Mq, or 50 gas 
particles. 

The number of galaxies within the virial radius of each 
halo at the present day is plotted against the halo mass in 



Fig. 10. The largest cluster contains nearly 30 large galaxies 
and has a virial radius of over 1.5ft Mpc. Its galaxy con- 
tent is consistent with those of 10 15 /i _1 Mq clusters which 
typically contain 30 - 100 L* galaxies. This is, in fact, the 
only large Abell-type cluster that formed in our simulations. 
All dark halos in the simulations that contain more than a 
few hundred particles have at least one galaxy within them. 
There are no empty halos of mass greater than 10 13 M© in 
either cosmology. 



4.2 Star formation and the age of galaxies 

Although, for simplicity, we have not included any prescrip- 
tion for star formation in our simulations, we can still de- 
rive some general conclusions regarding the properties of the 
stars expected to form. Stars form from gas that has cooled 
into galactic dark halos. The approximate time at which an 
individual gas particle is available for star formation is easy 
to estimate in the simulation. We take this to be halfway 
between the output time at which the particle first has an 
overdensity greater that 10 and a temperature smaller than 
12000 K, and the previous output time. These criteria have 
previously been shown to pick out solely baryonic material 
that has cooled into galactic objects (e.g. Kay et al. 2000). 
We assume that as soon as a particle becomes available for 
star formation, a mass of stars equal to the mass of the 
particle forms. 

By calculating a star formation time for each particle 
that ends up in galaxy, we can derive a mean, mass- weighted 
stellar age or formation time, tform, for each galaxy as 



JV 

N' 

i=l 



(1) 



where the sum is over all particles found within the final 
galaxy. This definition of formation time is based exclusively 
on the age of the stellar population of a galaxy and takes 
no account of when the galaxy was assembled. An alterna- 
tive definition of formation epoch is the time at which half 
the final galaxy mass was assembled into a single progeni- 
tor. As illustrated in Fig. 11, these two times can be very 
different. For the largest galaxy in the SCDM simulation vol- 
ume (which has a gas mass of 2.4 x 1O 12 M0 and for which 
we have enough resolution to follow many progenitors), the 
average formation time of its stars corresponds to redshift 
z ~ 2.5, but half of the stellar material is only assembled 
into a single distinct object at z ~ 0.7. The formation history 
of this galaxy bears little resemblance to the classic spheri- 
cal top-hot collapse model. Rather, the galaxy is assembled 
through repeated mergers of sub-units, many of which have 
a mass close to the resolution threshold of the simulation 
and which are originally spread over quite a large comoving 
volume. These fragments tend to line up with filaments of 
the dark matter distribution, giving rise to an anisotropic 
fragment distribution which is particularly visible at z = 3 
in Fig. 11. The mass distribution of the progenitors of this 
massive galaxy at various epochs is shown in Fig. 12. 

Fig. 13 shows the distribution of the mean stellar ages 
(as a fraction of the lookback time) for galaxies in the SCDM 
simulation. Contrary to the naive expectation for a hierar- 
chical scenario, there is a weak trend for the largest galaxies 
to be the oldest. This may be an artificial consequence of 
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Figure 10. The number of galaxies with mass above 6.4 X 10 10 Mq found within the virial radius of each dark matter halo. In the 
ACDM model, the largest halo contains 29 galaxies (24 in SCDM) and there are 6 clusters of 10 or more galaxies (7 in SCDM). The 
largest clusters have suffered significant depletion of their galaxy population due to ablation of galaxies moving through the cluster halo. 




Figure 12. The mass distribution of progenitors of the most mas- 
sive galaxy in the simulation. Each histogram refers to a different 
redshift as shown in the legend. 



the resolution limit of our model which is at a fixed mass. 
Models with more realistic feedback and no lower mass limit 
produce a weak trend in the opposite sense for galaxies in 
the field and no trend at all for galaxies in clusters (Kauff- 
mann 1996). 

An estimate of the star formation rate (SFR) in our sim- 
ulated galaxies can be derived from the gas mass that cools 
between successive output times. This crude measure (we 
only have 10 output times) is shown for the 20 largest galax- 




5x10' 
Object Mass (M Q ) 

Figure 13. The age distribution of mean stellar ages of galaxies 
in the SCDM simulation. The age is shown as a fraction of the 
lookback time (i.e. corresponds to the present and 1 to the 
Big Bang). The largest galaxies are the oldest because their stars 
formed preferentially in small fragments at at high redshift (see 
fig. 11). 



ies selected at the end of the SCDM simulation in Fig. 14. 
Clearly, the inferred aggregate SFR within each of these ob- 
jects is very high but, as Fig. 11 shows, these galaxies are 
broken up into many small precursor objects at early times 
each of which has a modest SFR. The aggregate SFR for 
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Figure 11. The positions of the progenitor fragments that make up the largest galaxy at z = in the SCDM cosmology. Objects 
containing more than 10 particles are shown as open squares. At each redshift, the circle indicates the virial radius of the largest 
progenitor of the dark matter halo hosting the largest galaxy, and the grey shading shows the underlying dark matter density. Each 
frame illustrates a comoving cubic region of side 12h _1 Mpc. Half the mass has cooled by z ~ 2.5, but it is not assembled into a single 

object until around z ~ 0.7. r-i i — i 
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Figure 14. The star formation rate for the 20 largest galaxies by 
mass selected at z = in the SCDM simulation as a function of 
cosmic time. 



Figure 15. The star formation rate for every 50th galaxy ordered 
by mass selected at z = in the SCDM simulation as a function 
of cosmic time. 



the most massive objects today is sharply peaked at z ~ 3. 
This contrasts with the star formation history of more typ- 
ical galaxies, illustrated in Fig. 15. This shows SFRs for 
every 50 th galaxy ordered by mass. These curves are more 
extended in time and the mean SFR for the average galaxy 
is much lower than the aggregate star formation rate of the 
largest objects. The SFRs generally decline towards low red- 
shift as galaxies consume the gas that is able to cool by the 
present day. The SFR integrated over all the galaxies in the 
simulation is shown in Fig. 4. 

By convolving the SFR in each galaxy with the Bruzual 
& Chariot (1993) stellar population synthesis model, we 
have calculated the if— band luminosity of each galaxy. We 
have assumed a Salpeter Initial Mass Function with a low 
mass cutoff of 0.1 Mq and solar metallicity. (This is higher 
than our assumed metallicity of 0.3 solar for the ambient gas, 
since the galactic material will have been enriched by repro- 
cessed metals.) As shown in Fig. 16, the model luminosity 
function in both models has a similar shape to the K— band 
luminosity function measured by Gardner et al. (1997). Fur- 
thermore, the model luminosity function can be made to 
agree quite well with the data by assuming that only a frac- 
tion 1/T of the stellar mass is luminous (and the rest is 
in brown dwarfs), with T = 2.8 for ACDM and T = 0.83 
for SCDM. The former is somewhat larger and the latter 
somewhat smaller than the values inferred from observa- 
tions (Cole et al. 2000 and references therein) because, as 
discussed earlier, we have used the same mass per gas par- 
ticle in the two simulations and, as a result, more gas cools 
globally in the ACDM model than in the SCDM case. Due 
to limited resolution, only the bright end of the luminosity 
function (L>L t ) is accessible in our simulations. 



5 CLUSTERING PROPERTIES OF THE 
GALAXIES 

In this section we examine both the two-point spatial clus- 
tering and the pairwise velocity dispersions of the galaxies 
as a function of separation and we compare them to obser- 
vational determinations. We also investigate the dependence 
of spatial clustering on galaxy mass. We test a virial mass 
estimator, that of Heisler, Tremaine, & Bahcall (1985, here- 
after HTB), on the galaxy samples contained within a few of 
the largest dark matter halos to see how well the true mass 
within the virial radius is recovered. 

Before proceeding with our analysis, we need to define 
a suitable galaxy sample. As discussed earlier, the galaxies 
are easily identified because of their very high gas density 
compared to the mean. The exact membership of the galaxy 
catalogues obtained using a group finder is relatively insen- 
sitive to the input parameters used to select them. In this 
section, galaxy catalogues for both cosmologies were gen- 
erated by applying a "friends-of-friends" group finder with 
a linking length of 0.02(1 + z) on a subset of gas particles 
with temperature less than 12000 K. All gas particles were 
counted for the purposes of defining the linking length itself. 
Catalogues were constructed at z~0, 1 and 3. They contain 
2263, 2360 and 1502 galaxies respectively for ACDM and 
2089, 1647 and 550 galaxies for SCDM, counting objects of 
8 particles or more. 

5.1 The galaxy autocorrelation function 

Figure 17 shows the redshift evolution of the two-point corre- 
lation function of galaxies in the ACDM and SCDM models. 
In both cases, the galaxy correlation function evolves little 
over the range z =0-3, in contrast to the mass autocorrela- 
tion functions also shown in the figure. The galaxy correla- 
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Figure 16. The K— band luminosity function. The left panel shows the ACDM model and the right panel the SCDM model. The open 
squares show our model results, while the filled squares are observational data from Gardner et al. (1997). Poisson error bars are shown 
in both cases. 



tion functions are essentially unbiased on the largest scales 
relative to the z = mass correlation function but show a 
distinct anti-bias on scales below 3/t -1 Mpc for ACDM and 
2 ft -1 Mpc for SCDM, and a positive bias on scales below 
~ 100 ft -1 kpc in both cases. Also plotted is Baugh's (1996) 
measured galaxy autocorrelation function determined by de- 
projecting the angular correlation function of galaxies in the 
APM catalogue. (The curve plotted assumes that cluster- 
ing is fixed in comoving coordinates.) On small scales, the 
galaxy correlation function in the simulations is encourag- 
ingly close to the real galaxy correlation function. On large 
scales, the SCDM model fails because it lacks sufficient large- 
scale power, but the ACDM model continues to do well. The 
departure of the ACDM correlation function of both galax- 
ies and dark matter below the APM result at the largest 
scales plotted is due to the small volume of the simulation. 
A better comparison between the APM data and the dark 
matter correlation function at large separations is shown in 
Fig. 5 of Jenkins et al. (1998). 

To examine the dependence of spatial clustering on 
galaxy mass, we subdivided the ACDM z — sample into 
low-, middle- and high-mass subsamples, each with approx- 
imately 750 galaxies. In terms of the gas particle num- 
ber, N P art, the three subsamples have 8 < N par t < 50, 
50 < Npart < 95 and N par t > 95 respectively, corresponding 
to if -band luminosities in the ranges (—20, —22.3), (—22.3, 
-22.8) and brighter than -22.8 . Figure 18 shows the 2- 
point correlation function for the three galaxy subsamples. 
The clustering strength shows a clear trend, with the most 
massive galaxies having the strongest clustering. The differ- 
ence between the medium- and low-mass samples is much 
less pronounced than the difference between the medium- 
and high-mass samples. 

To conclude, we find that the clustering properties of 
the galaxies in the simulations are markedly different from 



those of the dark matter. In particular, the galaxy correla- 
tion function evolves weakly with redshift and, coinciden- 
tally, appears unbiased on large-scales at the present epoch. 
At the present day, the correlation amplitude of the galaxies 
in both models is weaker than that of the dark matter at in- 
termediate scales around 1 Mpc and this is particularly 
pronounced for ACDM . The galaxy correlation function in 
the ACDM model is close to that determined observation- 
ally. Its shape is rather noisy because of the small sample, 
but it is much better fit by a power-law than the mass dis- 
tribution. The clustering strength of galaxies increases with 
galaxy luminosity, but this effect is only strong at very bright 
luminosities. 



5.2 The galaxy projected pairwise velocity 
dispersion. 

In this subsection, we compare the projected pairwise ve- 
locity dispersions of galaxies and dark matter. As is well- 
known, a reliable estimate of this statistic requires a rather 
large sample volume (Marzke et al. 1995, Mo, Jing & Boerner 
1997). The sample volumes that we analyse in this paper are 
probably too small to provide an accurate estimate of the 
global pairwise velocity dispersions, and we will therefore 
concentrate on the relative dispersions of galaxies and dark 
matter. Nonetheless, with this caveat in mind, we plot, for 
interest, the pairwise velocity dispersions determined from 
the LCRS catalogue by Jing, Mo & Boerner (1998). 

Fig. 19 shows the measured projected pairwise ve- 
locity dispersions of both dark matter and galaxies in 
the ACDM and SCDM models. The detailed definition 
of the projected pairwise dispersion is given in Jenk- 
ins et al. (1998). Because of our smaller simulation vol- 
umes, we have decreased the length-scale of the projection 
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Figure 17. Evolution of the galaxy two-point correlation function in the ACDM and SCDM models. Also plotted is the observed 
galaxy correlation function as determined by Baugh (1996) from deprojecting the angular correlation function of galaxies in the APM 
survey. (These data assume no evolution of the clustering with redshift, as is seen for the simulation galaxies). In both models, the 
galaxy correlation function in comoving coordinates changes little between redshift 3 and the present, in contrast to the mass correlation 
function which changes considerably. For ACDM the galaxies show a significant anti-bias on scales 0.2-3 h -1 Mpc and similarly for SCDM 
up to ~2h~ 1 Mpc . See text for more discussion. 



from ±25 h' 1 Mpc to ±10 h^ 1 Mpc and ±8 h' 1 Mpc for the 
ACDM and SCDM models respectively. We plot the same 
LCRS data points (taken from Jing et al. 1998), as plotted 
in Fig. 11 of Jenkins et al. (1998). 

There is a remarkable difference between the pairwise 
velocity correlations of galaxies and those of the dark mat- 
ter on scales below ~ 2 h' 1 Mpc in both cosmological mod- 
els. For ACDM , the galaxy pairwise velocities turn out to 
be not much higher than the observational data points. In 
fact, the dark matter pairwise dispersion in our 70 h' 1 Mpc 
ACDM simulation cube is around 100km s _1 higher than the 
value obtained by averaging over larger 239.5 h" 1 Mpc boxes 
in Jenkins et al. (1998; Fig. 11). One might infer from this 
that the global average for the projected galaxy pairwise 
dispersion is also overestimated. The difference between the 
behaviour of the galaxies and dark matter becomes smaller 
with increasing pair separation and is similar on the largest 
scales probed by these simulations. 

Why do the galaxies have a lower pairwise dispersion 
than the dark matter? To help answer this question we have 
constructed a "shadow galaxy" catalogue for ACDM by se- 
lecting the nearest dark matter particle to each galaxy. The 
shadow catalogue has, by construction, virtually an identi- 
cal 2-point correlation function to the galaxies. In Fig. 19 
we also plot, as a dotted line, the projected pairwise dis- 
persion for the shadow catalogue. This is almost identical to 
the projected pairwise dispersion of the galaxies. This shows 
that the strong difference between the galaxies and the dark 
matter arises mostly from the way in which galaxies popu- 
late dark matter halos rather than from a strong intrinsic 
bias in the velocities of galaxies relative to the dark matter 
particles at the same location. The fact that the dispersion 



of the shadow catalogue is everywhere higher than that of 
the galaxies indicates that some residual is present, but the 
closeness of the curves shows that this is not a large effect. 
Thus, the low galaxy velocity dispersion is a statistical ef- 
fect which may, in principle, be driven by differences in the 
ratio of the number of galaxies to total mass in high- and 
low-mass halos and also by differences in the way in which 
galaxies are distributed within halos compared to the dark 
matter. 

Our results for ACDM are very similar to those in Ben- 
son et al. (2000a), but differ in some respects from those of 
Kauffmann et al. (1999b). In both papers, semi-analytical 
modelling of galaxy formation was used to create synthetic 
galaxy catalogues from the same ACDM N-body simula- 
tion, but the detailed placement of galaxies in halos was 
different. Generally, in the semi-analytic models the effi- 
ciency of galaxy formation per unit mass of dark matter 
is a strong function of dark halo mass and peaks for halos 
with masses around 1O 12 M0. Benson et al. (2000a) found 
that by populating the dark matter halos unevenly as their 
semi-analytical model predicts, not only does one obtain a 
galaxy correlation function which matches the APM result 
of Baugh (1996) well, but also a comparable difference as 
found here in the projected pairwise dispersion of galaxies 
relative to the dark matter. Kauffmann et al. (1999b) find a 
very similar two-point correlation function, but rather little 
difference between the pairwise dispersions of galaxies and 
dark matter. As the analysis of Benson et al. (2000a) shows, 
this discrepancy can be traced back to relatively small varia- 
tions in the precise form of the halo occupation number pre- 
dicted in the two models, particularly for halos with mass 
> 1O 13 M0. These generate a large difference in the pairwise 
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Figure 19. A comparison of the projected two-point pairwise velocity dispersion of galaxies and dark matter for ACDM and SCDM. The 
data points are from Jing et al. (1998) based on the Las Campanas Redshift Survey. The galaxy-galaxy pairwise velocities are significantly 
different to those of the dark matter. The dotted line shows the projected two-point pairwise velocity dispersion of a catalogue made by 
selecting the nearest dark matter particle to each galaxy. The similarity of these dispersions and those of the galaxies suggests that the 
differences between galaxies and dark matter as a whole are mostly due to the way in which galaxies populate dark halos rather than to 
an intrinsic velocity bias of galaxies. 



dispersions but have a much weaker effect on other statistics 
such as the two-point correlation function. A detailed com- 
parison between our simulations and semi-analytical mod- 
elling techniques applied to dark matter only versions of our 
simulations is given in Benson et al. (2000c). 

In summary, we find significant differences in the pair- 
wise velocity dispersion of galaxies relative to the dark mat- 
ter on scales below 2 Mpc . The amplitude of this statis- 
tical effect appears sufficient to reconcile ACDM with cur- 
rent observations. The differences in the velocity dispersions 
of galaxies and dark matter arise because the galaxies pop- 
ulate dark halos with a variable efficiency which decreases 
for the highest mass halos. 



5.3 Virial theorem mass estimates of Galaxy 
clusters 

In this section we employ a commonly used virial estimator 
to see how well the virial masses can be determined for a 
few of the most massive dark halos using the projected po- 
sitions and velocities of the galaxies located within the halo. 
The estimator we have selected is one of those discussed by 
Heisler et al. (1985), specifically, 



Mvt = 



3nN J2> V h 
2G E^-l/Jt 



(2) 



where V r ,i is the line-of-sight velocity of the ith galaxy (in 
the observer frame where the cluster is at rest) and Rij is the 
projected distance between galaxies i and j. HTB show that 
this estimator is fairly robust. It does assume, however, that 
the galaxies trace the dark matter. If the galaxies are, for in- 



stance, more concentrated within a dark halo than the dark 
matter then the kinetic energy-like upper sum in equation 
(2) will be underestimated, whilst the potential energy-like 
sum will be overestimated. The result then is a virial mass 
estimate which will be lower than the true mass. Additional 
sources of uncertainty in applying the estimator to galaxy 
clusters are that the dark halos are not that well isolated 
from their surroundings and may not be in close virial equi- 
librium particularly if viewed at an epoch of major mass 
accretion. 

As a sample, we take the 9 most massive clusters from 
both ACDM and SCDM models. The average number of 
galaxies in our ACDM halos is N=12.4 and N=10.7 for 
SCDM. We project each cluster from a large number of dif- 
ferent directions to get the full distribution of the estimated 
virial mass. We then average the ratio of the estimated mass 
to the true mass for the cluster samples. Using log 10 (M) as 
the mass variable we find that the distribution of the es- 
timated mass about the true mass is very similar in both 
cosmologies but it is slightly more peaked in ACDM than 
SCDM. With only nine halos for each model, the conclu- 
sions that we can draw are not very strong. The distribu- 
tions have significant tails towards low masses, but the exact 
form is dependent on only one or two of the halos. The esti- 
mator shows a bias in log 10 (M) of 0.08±0.08 for ACDM and 
—0.02 ±0.2 for SCDM which is consistent with zero. A more 
reliable statistic which does not depend so much on the tail 
of the distributions is the range ±Alog 10 M about the true 
mass which includes 3/4 of the distribution. For ACDM this 
is about 0.23 and for SCDM 0.2. In Figure 20 we combine 
the results for ACDM and SCDM. For the combined sam- 
ple, 3/4 of the distribution is included within a range of 
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Figure 18. The ACDM galaxy correlation function for galaxies of 
different mass (or luminosity). The ACDM catalogue was sorted 
by mass and divided into three approximately equal sample sizes 
(see text for details). The most massive galaxies exhibit stronger 
clustering than the low-mass galaxies. 

±0.21. This spread of values is a little larger than the cor- 
responding spread found by HTB, who quote a range of 0.2 
for N=5, and 0.15 for N=10 but these are for realisations 
of virialised Mitchie models. This is not perhaps surprising 
given that the galaxies do not trace the mass exactly nor 
are the clusters in complete virial equilibrium. 

In Figure 21 we compare the average radial distribu- 
tion of the galaxies and mass within the virial radius for the 
9 halos from each cosmology analysed above. The galaxy 
radial distribution is shown with a solid line and the mass 
radial distribution with a dotted line. There are 112 galaxies 
contributing to the ACDM histogram and 97 to SCDM. All 
the histograms are normalised so as to have unit area under- 
neath. Broadly speaking the galaxies and mass have a similar 
distribution. A \ 2 analysis of the galaxies and mass distribu- 
tions yields \ 2 - 14 for 8 d.o.f. for both ACDM and SCDM. 
Taken individually this is unremarkable with a higher value 
occurring approximately 8% of the time. The joint x 2 nas a 
probability which is only marginally more significant at 3%. 
Without a larger sample it is difficult to draw any firmer con- 
clusion about the relative distributions except to note that 
several hundred galaxies are needed in order to distinguish 
the distributions. 

In summary, we find that applying the estimator, eqn. 2, 
to the galaxies in the largest dark matter halos taken from 
either of our two simulations, one can infer the masses of the 
dark halos in a relatively unbiased way, albeit with a signif- 
icant uncertainty for a single determination of a particular 
halo. Consistent with this conclusion, the radial distribu- 
tions of galaxies and mass are similar. These mass estimates, 
of course, neglect projection effects in the identification of 
cluster galaxies which can add significantly to the uncer- 
tainty (eg van Haarlem, Frenk & White 1997). It is interest- 



Figure 20. Combined distribution function for the virial mass 
estimator (eqn. 2) applied to the galaxies within the largest 9 
dark matter halos in our ACDM and SCDM simulations. The 
x-axis shows the difference between the estimator and the true 
mass within the virial radius. The estimator is computed over a 
large number of randomly chosen directions about each halo. The 
average mass of the combined samples is consistent with the true 
virial mass; three-quarters of the distribution lies within ±0.21. 
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Figure 21. The average radial distribution of galaxies (solid 
lines) and mass (dotted lines) within the virial radius of the 9 
largest halos from ACDM and SCDM simulations. All the his- 
tograms are normalised to have unit area underneath. The distri- 
butions arc similar. 
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ing to note that for clusters more massive that ~ 1.3 x 10 
Mq, the mean baryon mass fraction in the form of galaxies 
is 17% and 8.4% for ACDM and SCDM respectively. These 
values differ from the global fractions which are 12% for 
ACDM and 7% for SCDM. Thus, the mass-to-light ratios 
of clusters cannot, in general, be used to estimate S7o reli- 
ably without accurately accounting for the difference in the 
efficiency of galaxy formation in clusters compared to the 
universe as a whole. 



6 DISCUSSION AND CONCLUSIONS 

Our simulations have many limitations and should be re- 
garded primarily as illustrative of general trends likely to 
characterize galaxy formation in hierarchical CDM models. 
Most of these limitations stem from the limited resolution 
that is attainable with current computing resources if one 
wishes to simulate a representative volume of the universe. 
Our models do not include a self-consistent treatment of star 
formation and feedback. Instead, the ability of gas to cool in 
small objects is limited by the resolution of the simulation. 
This is clearly a crude approximation. The one virtue of the 
strategy we have adopted is that the resolution limit can be 
adjusted to ensure that approximately the right amount of 
gas cools by the end of the simulation. We have done this 
for the SCDM simulation and kept the same resolution for 
the ACDM simulation with the result that more gas cooled 
in the latter case than is observed in the form of stars in 
the real universe. Clearly, further progress will require more 
realistic, resolution-independent and physically-based treat- 
ments of star formation and feedback. 

Another limitation of our simulations, stemming also 
from limited resolution, in this case spatial resolution, is the 
inadequate modelling of the internal structure of galaxies. 
It seems clear that a gravitational softening of 10/t -1 kpc is 
much too large to prevent the artificial disruption of galax- 
ies in rich clusters. Our strategy of retaining gas that has 
cooled into galaxies as a gaseous component rather than 
turning it into stars as is sometimes done, partially com- 
pensates for this, but it is clear that simulations with much 
higher spatial resolution are required. Unlike the problem of 
star formation and feedback, this is a limitation that can, in 
principle, be overcome by increased computing power. For 
example, a simulation similar to ours but with a gravita- 
tional softening of only 2h~ l kpc would require ~ 100,000 
timesteps at least for the particles in the densest regions. 

In spite of these limitations, our simulations demon- 
strate the potential of direct simulation for realistic mod- 
elling of galaxy formation. This study is the first to follow 
the formation and evolution of a large number of galaxies in 
a volume, a cube of side 100 Mpc, which may be regarded as 
representative. Over 2000 galaxies brighter than L* formed 
in each of our simulations. Both in this paper and elsewhere 
(Pearce et al. 1999) we have shown that the galaxies that 
form in the simulations have a spatial distribution that is 
consistent with observations not only locally, but also at 
high redshift. 

Here, we have used our simulations to illustrate the hi- 
erarchical build up of galaxies and to give an indication of 
the rate at which gas is expected to cool into fragments 
which subsequently merge to form larger galaxies. The ag- 



gregate star formation rate in the fragments that end up in a 
massive galaxy today adds up to a large value and is sharply 
peaked at high redshift. By contrast, the star formation rate 
of a typical galaxy today is more modest and spread out in 
time. The mean star formation rate in the simulation as a 
whole rises at early times, reaches a broad maximum be- 
tween z — 2 and z = 1 and declines towards the present, 
in a manner reminiscent of the data of Steidel et al. (1999). 
Although this behaviour is, to a large extent, determined by 
resolution effects, the fact that it resembles the data suggests 
that the bulk properties of the model galaxy population to- 
day are not too unrealistic. In the ACDM simulation, the 
first substantial galaxies form at z ~ 5. The galaxy popu- 
lation builds up rapidly until z = 1 and there is a marked 
decline in the rate of change of the galaxy mass function af- 
ter z — 0.5. In the SCDM simulation there is more evolution 
at recent times. All these trends have been previously em- 
phasized in semi-analytic models of galaxy formation (e.g. 
Cole et al. 1994, Kauffmann 1996). At the present day both 
simulations match the bright end of the observed K-band 
luminosity function for suitable values of the mass-to-light 
ratio. 

The galaxy autocorrelation functions evolve little with 
redshift over the range z — — 3. The galaxy correlation 
function agrees closely with that of the dark matter at sep- 
arations of 10/i -1 Mpc at z = but differs significantly on 
smaller scales. For both, the ACDM and SCDM models, the 
galaxy correlation function appears closer to a power law 
than the mass correlation function. The amplitude of the 
galaxy correlation function increases with the mass of the 
galaxies. The projected pairwise velocity distribution of the 
galaxies is significantly lower, particularly in ACDM , than 
that of the dark matter. 

Within the most massive dark matter halos the galax- 
ies trace the dark matter faithfully. A virial mass estimator 
applied to the clusters correctly infers the appropriate dark 
matter mass albeit with a large dispersion for a single de- 
termination. However, galaxy formation is more efficient in 
such clusters than in the simulation as a whole, suggesting 
that the cosmic density parameter cannot safely be inferred 
from the mass-to-light ratios of clusters and the mean lumi- 
nosity density of the Universe. 
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